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We present a lattice QCD calculation of the A (1232) matrix elements of the axial- vector and 

pseudoscalar currents. The decomposition of these matrix elements into the appropriate Lorentz 

^1^ , invariant form factors is carried out and the techniques to calculate the form factors are developed 

I ■ and tested using quenched configurations. Results are obtained for 2 + 1 domain wall fermions 

\^0 ' and within a hybrid scheme with domain wall valence and staggered sea quarks. Two Goldberger- 

^h , Treiman type relations connecting the axial to the pseudoscalar effective couplings are derived. 

These and further relations based on the pion-pole dominance hypothesis are examined using the 

lattice QCD results, finding support for their validity. Utilizing lattice QCD results on the axial 

charges of the nucleon and the A, as well as the nucleon-to-A transition coupling constant, we 

perform a combined chiral fit to all three quantities and study their pion mass dependence as the 

chiral limit is approached. 
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PACS numbers: 

I. INTRODUCTION 

Great progress has been made in lattice QCD studies of hadron spectroscopy and structure and lattice QCD results 
are beginning to provide input to phenomenology and experiment. Simulations with dynamical quarks near and at 
the physical pion mass ll|-y| have been shown to produce the observed low-lying hadron spectrum [2, S |6[ and tt^ — 7r+ 
scattering lengths have been calculated to good accuracy j7Hl(j|. 

Whereas producing experimentally measured quantities from first principles provides a powerful validation of the 
lattice QCD methodology, calculating quantities that are difficult to extract or have impact in probing physics beyond 
the standard model is a much more challenging prospective. Studying the structure of the A resonance is an example 
of the input lattice QCD can provide to phenomenology that cannot be directly extracted from experiments. This 
is because the A decays strongly with a lifetime of ^ lO"^"^ seconds [ll|j [l2l and resists experimental probing. 
Measurements of the A"*" magnetic moment exist albeit with a large experimental uncertainty. The A, having width 
r ~ 118 MeV and lying close to the ttN threshold, plays an important role in chiral expansions. In heavy baryon chiral 
perturbation theory it has been included as an explicit degree of freedom [13l4l6| , where it is argued that it improves 
chiral expansions applied in the description of lattice QCD results such as the nucleon axial charge [Sa]. Chiral 
langragians with A degrees of freedom involve additional coupling constants that are difficult to measure. Therefore, 
one either treats them as free parameters to be fitted along other parameters using lattice QCD results [i3,li3| and data 
extracted from partial- wave analysis of scattering measurements [15l . Ilq or estimates them based on phenomenology 



and symmetries. For example, one can relate the nucleon axial charge qa, which is well measured, to the A axial 
charge, in the large- iVc limit ^^ or using SU{4:) symmetry [1^. The Goldberger-Treiman (GT) relation is then used 
to get the effective ttAA coupling. Another framework to extract the ttAA coupling is via sum rules |2C| . 

Lattice QCD provides a nice framework to study the A properties and calculate the A coupling constants. In some 
of our recent work we developed the formalism to study the N~ A transition form- factors within lattice QCD [2lj, ^^ , 
as well as the A electromagnetic form- factors [23 • The quadrupole electromagnetic form factor, extracted for the first 



time, provided input for the deformation of the A showing that in the infinite momentum frame the A is prolate [24 1. 

In this work, we present a detailed study of the axial- vector and pseudoscalar form factors of the A. The theoretical 
framework and a subset of the results were given in Ref. [25| . Here we discuss in detail the lattice techniques developed 
and utilized for the extraction of these form factors. In addition, we present an extended analysis of the momentum 
dependence of all the form factors using an additional ensemble of dynamical domain wall fcrmions. We also include 
a study of the pion-pole dominance predictions and compare them to our lattice QCD results. 

The outline of the paper is as follows: In Section II we present the decomposition of the A matrix elements of the 
axial-vector and pseudoscalar currents. In Section III we explain our lattice techniques and discuss the ensembles 
utilized for the calculation. In Section IV we present the lattice results on all form factors and examine several 
relations among them and their phenomenological consequences. In Section V we perform a combined chiral fit using 
our results on the nucleon axial charge qa [26[, the A axial charge Gaa calculated in this work and the dominant 
axial N — A transition form factor, C^, calculated in previous work on the same sets of lattices ^^. Finally, in 
Section VI we give a summary and conclusions. Technical details and our values on the form factors are presented in 
the Appendices. 

II. THE AXIAL AND PSEUDOSCALAR MATRIX ELEMENT OF THE A 

Lorentz invariance and spin-parity rules determine the decomposition of the A+ matrix element of the isovector 
axial- vector current in terms of four invariant functions of the momentum transfer squared, q^ = {pf — PiY". 

{A+{pf,Sf)\A^mA+{p,,s.,))=u^{pf,sf) [O^'^Y^ u^[p.,s,) 
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where Pi{si) denotes the initial momentum (spin) of the A and Pf{sf) the final momentum (spin). The flavor- isovector 
axial- vector current operator is defined as 

A^(x)=V^(x)7'^75y^(a:) (2) 

where t^ denotes the Pauli matrix acting in flavor space and ip{x) is the isospin quark doublet. The four axial form 
factors, gi,g3,hi and h^ as defined in Eq. ([T]) are grouped into the familiar structure of the nucleon axial-vector 
vertex. 

In the description of spin-3/2 energy-momentum cigenstatcs, classical solutions of the Rarita-Schwinger equation 
play a central role. Each component of a vector-spinor u^r, with a a Lorentz four- vector index solves the free Dirac 
equation 

y-MA]u^ip,s) = 0. (3) 

Implementing additionally the constraint equations, 

p^M^(p,s) = and 7"«^(p,s) = 0, (4) 

the unphysical components are eliminated and the remaining eight degrees of freedom describe a spin-3/2 (anti- 
)particlc. Rarita-Schwinger spinors satisfy the spin sum relation: 

3/2 

Aar = X] ^a{P^s)u^{p,s) 
s=-3/2 
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where the normalization u u'^ = — 1 is assumed. 

The zero momentum transfer Hmit of ([1]) defines the axial charge, Gaa of the A multiplet. Ref 231 normalizes the 
axial charge via 



(A++|4|A+^ 



{^-\Al\A-) 



GaaM, 



where A^^ encodes the spin structure of the forward matrix element 

Following the above normalization we establish via Eq. ([1]); 

Gaa = -351 (0). 



(6) 



(7) 
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We note that while the Lorentz decomposition of the axial current is naturally expressed via gi,g3,hi and h^ as 
in Eq.([T]), a decomposition in terms of multipoles is possible as for example in the case of the A electromagnetic 
transition |23| . Such a representation is more easily expressed in the Breit frame. This decomposition is performed 
in Appendix A in terms of four multipoles, ii, L3, Ei and E^, and their relation to the form factors gi, gs, hi and ft.3 
is given. 

The A+ matrix element of the pseudoscalar density operator 



P{x) =-0(x)75y0(a;) 

is decomposed in terms of two Lorentz invariant form factors, denoted by g{q^) and h(q^): 

{A+ipf,sj)\PmA+{p,,s.,)) ^u^{pf,sf) [O^Y^uHp^,s.), 



(9) 
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While g{q^) and h{q'^) are the directly computable form factors from the three-point pseudoscalar correlator, they 
can be related to the phenomenologically more interesting pion-A vertex using the partially conserved axial current 
hypothesis (PCAC). Using PCAC on the hadronic level one can write 



df^Al = UtuIt:^ , 



(11) 



with 7r° denoting the isotriplet pion field operator. In the SU(2) symmetric limit of QCD with m, denoting the 
up/down mass, the pseudo-scalar density is related to the divergence of the axial-vector current through the axial 
Ward-Takahashi identity (AWI) 



d^A''^ = 2m,P- - Umln" 



(12) 



with operators now defined as quark bilinears. Using the relations of Eqs. pT|) and ([T^ we identify the physically 
relevant pion-A-A from factor G^aa('?^), which at g^ = gives the ttAA coupling, as well as a second form factor 
H-KAA{q'^)i by rewriting the pseudoscalar matrix element as 



2m,(A+(p/,5;)|P(0)|A+(p„s,)) 



Um^ 



xut 



where we effectively make the identification 
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(14) 
(15) 



At zero momentum transfer q^ = Q only G-nAA can be extracted. This coupling is analogous to the known tt — N 
pseudoscalar coupling constant Gt^nn defined for the nucleon. For the discussion presented in the next section it is 



useful to recall the definition of the corresponding quantities in the nucleon sector |2l| . For the matrix elements of 
the axial-vector current we have 



{N{pf,Sf)\Al\N{p,,s,))^i^UN 



GAiq'Ha, + I^G,{q') 



UN 



(16) 



and for the pseudoscalar density 

2mg{N{pf,Sf)\P''\N{p,,s,)) 



{^ - ml) 



X UN [GT,NN{q^)\ i-t^UM 



(17) 



Note that we have dropped for simplicity an overall kinematical factor arising from the normalization of lattice states, 
since it is of no relevance for our discussion here. 

A. Goldberger-Treiman Relations 

In this section we apply PCAC to derive GT relations for the A. We recall that PCAC has been shown to apply 
satisfactorily in the nucleon case leading to the Goldberger-Treiman (GT) relation. This can be derived from Eqs. P^ 
and (|17p related by AWI and taking g^ = to obtain G-^nn in terms of the nucleon axial charge via the relation 



UGT^NNiO) = rriNGAiO)- 



(18) 



Assuming G^^nn varies smoothly with q^ so that G',rArjv(0) ^ G-nMNifriD = g-^NN then the GT relates the physical 
coupling constant Qt^nn with the nucleon axial charge gA- At the chiral limit, using d^A^^ ~ one derives that 

Gp{q^) = — -t^Ga('?^). Therefore, g-^NN measures the chiral symmetry breaking. PCAC dictates that the form 

factor Gp{q^) has a pion pole given by Gp{q'^) = "i_l Gj^NNiq'^)- The validity of the GT relation and the momentum 

dependence of Gp{q^) in the nucleon case has been studied in Ref. [21|. Similarly, a non-diagonal GT relation, 
applicable to the axial iV-to-A transition is formulated and relates the axial A^A coupling ca to the irNA effective 
coupling. Lattice calculations examined the validity of the non-diagonal GT relation using the same ensembles as in 
this work [22 1. 

One can similarly derive GT relations for the A by taking the matrix elements of the AWI with A states, 
(A|9^^^|A) = 2TOg(A|P|A). Taking the dot-product of g^ with the matrix element of the axial-vector current 
given in Eq. ([T]) we obtain 



rriA 



q"q'^ 



{ml - g2) 



q^q^ 



4iW A 



(19) 



where r = — g'^/(2A/A)^. By considering ct ^ t in Eq. (|19p we derive the relation 



TTiA {hi — rh^ 



UmlH^AA{q'^ 
ml — (p 



(20) 



which implies that 



A/a (.91 - rgj) 



UmlG^AAJq'^) 
m,1 — q^ 



(21) 



One possible linear combination of Eqs. ((20)) and ([2T|) can be obtained by taking the dot product of Eq. (|T9|) with 
qr leading to 



f Tfi 

Ma [{gi - rgs) - T{hi - T/13)] = —^ — ^ [G^tAA - tHt^aa] 

m- - q^ 



(22) 



which can be considered as a generalized GT-type relation connecting all the six form factors. At g^ = and assuming 
all terms in Eq. ([^^ are finite we obtain 



/7rG7rAA(0) = mA5l(0). 



(23) 



If GttAA is a continuous slow varying function of q^ as q^ — > then G-^i\i:^(vr?^') ^ GttAaCO) and wc thus derive a GT 
relation for the A analogous to the one for the nucleon case. 

Using Eq. [201 and setting q^ = we obtain a second GT relation 

UH^AA{0) = mAhi{0). (24) 



If one invokes pion-pole dominance and noting that gi and GttAA are both finite at the origin, it follows from Eqs. ([21 
[^n|) that as q^ — > mj g^ and h^ must have a pole at q^ = m^. We thus arrive at the relations 



U ^ _ 4/^AfA 

A'/a m^ - q^ 



51 = TJ-GttAA , 53 = — 2"^-^G^AA ■ (25) 



and 



U „ , 4/.AfA 

7-;— ^ttAA , "3 = — 5 5 

A'Ya m^ - q 



hi = -j^Ht,aa , /^3 = ^-^f^-^H-nAA ■ (26) 



It is thus interesting to note how the spin-3/2 nature of the A state combined with PCAC leads to a pair of Goldberger- 
Treiman relations, given by Eqs. ([23| and ([24]) . Let us examine further these relations at the chiral limit. From Eq. ([22]) 
we find that 

hi-Ths^ , (27) 

r 

which means that in the limit q^ — ;> the leadin behavior of hi ~ l/g^, /is ~ l/(q^)^ via Eq. ([20]) and Ht^aa ^ l/g^ 
via Eq. ([26]) . Therefore, the second GT-type relation given in Eq. ([24]) cannot be extrapolated to physical pion mass 
since the assumption that hi and Ht^aa are slowly varying functions of q^ no longer holds. However, since they both 
display a pion-pole behaviour one can factor it out on both sides and thus the ratio /ii/i?7rAA can be extrapolated to 
the physical pion. In this sense, this constitutes a second GT relation. 

III. LATTICE EVALUATION 

A. Euclidean Correlators and Form Factors 

Standard techniques are employed on the Euclidean space-time lattice for the evaluation of hadronic form factors. 
The following two-point and three-point functions are required: 

G,.(^^p,</) = ^e-"/-p^r:;,„(x.a(t/,x/)x..'(o,o)) 

G^^,(^^q,i;t/) = '£e+^''-'^K,^{xao.itf,Xf)A^{t,x)Xr^,iO,6)) 

X,Xf 

Glf{T-',q,t-tf) = X]e+"-'"r:;,„(x.c.(t/,x/)P(i,f)x™.(0,0)), (28) 

X,Xf 

where P{t,x) and Af^(t,x) are the lattice pseudoscalar or axial current insertions, and x is the standard lattice 
interpolating field with overlap with the A^ quantum numbers: 

X^^ix) = ^£-"^[2 (u"^(x)G7.d^(x)) <(x) + (u'^^(x)G7.u^(x)) d^(x)] . (29) 

The overlap of x with the spin-3/2 A+ is 

{n\xao.mAip,s))^Zut{p,s), {Aip,s)\xao.iOm)^Z*utip,s). (30) 

We will use the following F matrices, which project onto positive parity for zero momentum, for our calculation 

r^ = ^(l + 7'), r'-^ = ^(l+7^)757fc, fc = l,2,3. (31) 



The Fourier transforms in (|28p enforce a static A sink at the timeshce tf and a momentum transfer q= —p injected 

via the operator insertion at an intermediate timeslice t. 

We insert into these correlators complete sets of hadronic energy momentum cigcnstates: 
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A/„ 



VE, 



n(p) 
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(32) 



where with ^ we denote collectively all quantum numbers including spin. For large Euclidean time separations t and 
tf — t the ground state propagation dominates the correlator: 
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+excited states 
The Wick-rotatcd axial and pseudoscalar operators take the form 



(33) 
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(35) 



with the Euclidean four-momentum transfer Q^ — (j(A/a — i*^A(p)): ^q)- The Rarita-Schwinger spin-sum relation 
becomes 



A 



2Ma 



lair , '2'PaPr .Pair - Prla 

I- 



Mil 
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(36) 
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where all the 7 matrices are in Euclidean space: 70 = 74 and 7^^ 
Forming an appropriate ratio of the 3-point to the 2-point correlator serves to cancel out the unknown Z-factors and 
leading time-dependence. A particular product of 2-point correlators which minimizes the denominator noise-level is 
utilized as it contains smaller time-extents. The proposed ratios are: 



"^ ' Gfefc(r4,o,i/)VGfcfe(r4,o,i/-i)Gfefe(r4,p*,i)Gfcfc(r4,K,</) 



and 



K^T^Q^t) 
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(37) 



(38) 



for the axial and pseudoscalar vertices. Summation over A: = 1,2,3 is implicit in the 2-point correlators. At large 
Euclidean time separations tf — t and t these ratios become time- independent (plateau region) . 



jX (j^i> r^ +^ ^ nviX _ ^^^ 



i?^(^),(^^Q,t)^G^ 



CT(/i)T 



r''A..,(0)Of(^),A.,.(p) 



(39) 



where X stands for the axial {A^) or pseudoscalar (P) current. It is easy to show that the 2-point correlators are 
dominated by 



o.(^^,<,H.^— ^^1^(1.4 



(40) 



and therefore the constant C is determined as 



C=\^ 



2i?A(pi) 



M, 



A 



Ml 



Mi 



Mi 



(41) 



There are at most 256 available combinations of the Dirac and Lorentz indices in Equation (|39)) , each one expressed 
as a linear combination of the axial (pseudoscalar) form factors times kinematical tensor coefficients. Since we are 
interested in the momentum dependence of the matrix elements, evaluation of the 3-point correlators is required for 
a large set of transition momenta q for both A^^ and P operators. In order to perform this economically we utilize 
the sequential inversion through the sink technique [23[ by fixing the sink timeslice tf and performing a backward 
sequential inversion through the sink. The sequential vector is coupled with a forward quark propagator and the Fourier 
transformed insertion operator at all intermediate time-slices < t < i/ at a small computational cost, obtaining thus 
the full momentum dependence of the amplitude. A drawback in this approach is the fact that the quantum numbers 
of the source and sink interpolators -which correspond to the Lorentz indices a, t and V^ are now fixed per sequential 
inversion. Within the space of 64 available 3-point correlators corresponding to choices of a, t and v we perform an 
optimization by forming appropriate linear combinations such as the degree of rotational symmetry of the summed 
correlator is maximal and consequently all transition momentum vectors q that correspond to a fixed virtuality q^ 
will contribute to the form factor measurement in a rotationally symmetric fashion. This optimization technique has 
proved extremely useful in obtainin g hi gh accuracy results in the Nuclcon elastic, Nucleon-to-A electromagnetic [22l |. 
axial and pseudoscalar transitions 2l[ as well as the A electromagnetic form factors [23]. We evaluate the Dirac 
traces in Eq. ([55]) using symbolic software such as form [2g| and Mathematica. 
We construct the following two optimal linear combinations, which we refer to as Type-I and Type-II. 



3 3 



Type - I : nl^{q)^Y.Y. '^-^r [VA,,, {Q)0^;^^,AMp) 
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3 

Type - II : UJ/^iq) ^ ^ T^.tr [r^A,,, iQ)OXr'^r'r{p) 



with the matrix T: 



Tar 



a.T — 1 



1 -1 

1 
-1 



(42) 
(43) 

(44) 



Detailed expressions for the decomposition of the above combinations to the four axial form factors are provided 
in Appendix B. The above types are in addition utilized for the extraction of the two pseudoscalar couplings (Ap- 
pendix B). A large number (0(10'^)) correlators of axial (j4^,/i = 1,2,3,4) and pseudoscalar (P) insertion momenta 
q are combined for momentum transfers ranging up to ~ 3 GeV^ per ensemble. We stress that only two sequential 
inversions through sink -one for each Type above- are required in order to disentangle completely all six form factors 
from the relevant 3-point functions. 

Correlators corresponding to a fixed momentum transfer g^ are analyzed simultaneously in an overconstrained 
system analysis for the extraction of the form factors. Typically 0(20 — 50) plateau averages for the optimal ratios 
given in Eq. ([55]) will contribute to the determination of the form factor for each Q^ value. A global x^-minimization 
amounts technically to the singular value decomposition of an TV x M over-complete linear system, with M unknowns 
(4 for the axial or 2 for the pseudoscalar) and N input data (the 0(20-50) plateau averages). Further details on 
this kind of analysis can be found in Ref [29| . Jackknife estimates are utilized for all levels of variance extraction on 
observables. 



B. Ensembles and Parameters 



In Table U we summarize the parameters and number of configurations for the ensembles used in this work. As can be 
seen, three sets are employed. These are the same as the ones we used previously for the study of the nuclcon axial 
form factors as well as the nucleon-to-A axial transition form factors. Therefore, these ensembles provide a complete 
calculation of the nucleon/A sector, allowing a direct extraction of low energy couplings from a combined fit. 

The gauge configurations used in the analysis include a set of quenched configurations on a 32'^ x 64, at /3 = 6.0, 
corresponding to a lattice spacing a = 0.092 fm with pion masses 560 MeV, 490 MeV and 411 MeV. The low statis- 
tical noise makes this ensemble appropriate for checking our lattice methodology and some of the phenomenological 



V Stat, m^ (GeV) tun (GeV) niA (GeV) 

Quenched Wilson fermions 

l3 = 6.0, a-^ = 2.14(6) GeV 



32^ X 64 


200 


0.563(4) 


1.267(11) 


32'^ X 64 


200 


0.490(4) 


1.190(13) 


32^ X 64 


200 


0.411(4) 


1.109(13) 



1.470(15) 0.1554 
1.425(16) 0.1558 
1.382(19) 0.1562 



Mixed action, a"' = 1.58(3) GeV 

Asqtad (am„.d/3 = 0.02/0.05), DWF (am„,d = 0.0313) 

20^ X 64 264 0.498(3) 1.261(17) 1.589(35) 

Asqtad (am„.d/3 = 0.01/0.05), DWF (am„,d = 0.0138) 

Wx 64 550 0.353(2) 1.191(19) 1.533(27) 

Domain Wall Fermions (DWF) 

m„,d/ms = 0.004/0.03, a~^ = 2.34(3) GeV 

'Wx 64 1428 0.297(5) 1.27(9) 1.455(17) 

TABLE I: Ensembles and parameters used in this work. We give in the first column the lattice size, in the second the statistics, 
in the third, fourth and fifth the pion, nucleon and A mass in GeV respectively. We did not do a full form-factor analysis 
on the 20^ x 64 mixed-action ensemble. Rather we merely determined the axial matrix element at g^ = (a much cheaper 
computation) for our axial charge chiral fits. 

relations. We apply Gaussian smearing at the source and sink in order to minimize the excited state contamination on 
the baryon correlators. The parameters a = 4.0 and n = 50 have been tuned to provide optimal overlap to a nucleon 
state l2l|. The source-sink separation is set at AT = 12a = 1.1 fm. In our previous studies involving the A such a 
time separation was found sufficient for ground state dominance. The second set consists of two ensembles that use 
two degenerate light and one strange {Nf = 2 -|- 1) Asqtad-improved dynamical staggered fcrmions generated by the 
MILC collaboration [S^. The strange quark mass is fixed to its physical value, the lattice spacing is set to 0.124 fm 
and the lowest pion mass is 353 MeV. Our calculation employs Domain Wall (DW) valence quarks with light quark 
mass tuned so as the pion mass matches the lowest pion mass obtained using staggered fermions. The extent of the 
fifth dimension of the domain wall action is set to L5 = 16a, which was demonstrated to provide minimal violations 
to the chiral symmetry properties of the domain wall fermion (DWF) operator. The source-sink separation is set to 
AT = 8a = 1.0 fm and Gaussian smearing is applied at the source and sink with APE smearing on gauge links that 
enter the smearing function applied on the interpolating fields. The parameters are given in Ref. [2l[. Finally, the 
third set is an Nf = 2-1-1 ensemble of DWF generated by the RBC-UKQCD collaborations [3l| with a lattice spacing 
a = 0.084 fm and the physical volume of (2.7 fm)'^ and pion mass of 0.297 MeV. The extent of the fifth dimension is 
also here T5 = 16a. It turns out that the residual quark mass introduced via the chiral symmetry breaking effects 
is arrives = 0.000665(3), or 17% of the bare quark mass. The smearing parameters for the interpolating fields are 



given in Ref. [22|. The sink-source time separation is set at AT = 12a = 1.01 fm. In order to increase the statistics 
at this lowest pion mass we use the coherent sink technique, employed in our study of the nucleon to A transition 
using the same ensemble ^^. The four quark sources are placed at timcslice ti — {i — 1)16, i = I,-- • ,4 for each 
configuration. Four forward propagators must be computed - each with a source at one of the time-slices. The A 
sinks are constructed at all four equally spaced time-slices i/(j) ~ ti + 12 and one sequential inversion is performed 
in order to construct the three-point correlator. Gauge invariance ensures that combining the sequential vector with 
each one of the forward quark propagators generated at ti, projects the appropriate A matrix element between ti 
and ti + 12 as the other cross-terms will average to zero. It has been shown in Ref. [22| that while statistics is thus 
multiplied by 4, the noise level is not raised above what expected from the four completely independent correlators 
which participate in the coherent sink. This means that we can reduce the error by a factor of two at the cost of one 
sequential inversion. Therefore the 1428 statistics given in the table correspond to 357 coherent sequential inversions 
per each type of combinations (see Eq. (|43)) ). 

IV. RESULTS 

A. Axial-vector and Pseudoscalar Form Factors 

In this section we present results on the A axial-vector and pseudoscalar form factors from the ensembles utilized in 
this work. The axial current is renormalized multiplicatively in all ensembles. Values for the renormalization constant 
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FIG. 1: Lattice QCD results for tlie gi axial form-factor. 
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FIG. 2: Lattice QCD results for the 33 axial form-factor. 



Za are provided in Table HIl 

In Figs. [U m [3] and |4] we show the results for the four axial form factors, gi, 33, hi and /13, respectively. All the 
results on these form factors are provided in Appendix C. The form factor gi is the dominant axial- vector form factor 
and the only one that can be extracted directly from the matrix element at Q^ = 0, determining the axial charge of 
the A. Based on PC AC and pion pole dominance we expect gi to be a smooth function of Q^, whereas hi and (73 
to have a pion-pole and /13 a double pion-pole behavior. Given that gi and hi are multiplied by Q^, whereas h^ is 
multiplied by Q* it is increasingly more difficult to resolve these form factors via the simultaneous overconstrained 
analysis of the measured matrix element of the axial-vector current, especially at small Q^ ~ a fact that is clearly 
reflected on the statistical error of the form factors shown in the figures. The results from the quenched ensemble, 
although based on the analysis of 200 configurations, have the lowest statistical noise and this is the primary reason 
for using them in this first calculation of the form factors. The statistical noise is more severe for the DWF ensemble 
at ttItt = 297 MeV for which results on /13 are too noisy to be useful and are omitted from plots. We do, however, 
include these numbers in the tables in the Appendix C for completeness. 

Figs. [5] and [6] show the pseudoscalar form factors g and h, respectively, as defined in Eq. (fTO)) . where the pion 
pole is explicitly written. The numerical values of these form factors are provided in Appendix C. As confirmed 
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FIG. 4: Lattice QCD results for the /13 axial form-factor. 



by the numerical results, g is the dominant pscudoscalar form factor showing a pion-pole dependence, whereas the 
subdominant form factor h shows a stronger Q^-depcndence consistent with a double pion-pole. In section 11 we 
already defined the physically relevant pion-A coupling G'ttAaI'ti^) factoring out the pion-pole and fixing coefficients 
via PCAC through Eq. [TH G^aa(9^) has a finite value at the origin, as can be seen in Fig.[7]where numerical results 
arc depicted. This value in fact defines the traditional strong coupling ^ttAA of the pion to the A state via 



SttAA = G.„AAiml) . 



(45) 



The secondary momentum-dependent coupling, H^rAAiQ^), is plotted in Fig. [51 The numerical results are consisted 
with a pion-pole divergence at small Q^ as expected from the analysis given in the previous section. The statistical 
error on this coupling is larger in particular at small Q^ since in the combined analysis the pscudoscalar matrix 
element is multiplied by a factor of Q"^ . 

Notice that the extraction of Gt^aa and Ht^aa from Eqs. (fT4)) and (fTS)) requires knowledge of the light quark mass 



rriq and the pion decay constant, /tt, on each of the ensembles. Calculation of /^ requires the two-point functions of 
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FIG. 5: Lattice QCD results for the g pseudoscalar form-factor. 
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FIG. 6: Lattice QCD results for the h pseudoscalar form-factor. 



the axial-vector current A'^ with local-smeared (LS) and smeared-smcarcd (SS) quark sources, 



C 



LS 



it)=Y,{^\T (A3(x,t)il(O,0)) \Q 



(46) 



(and similarly for Cgg), where A^{x,t) denotes the local operator and ^Kx, i) the smeared operator. The pion decay 
constant /tt is obtained from the pion-to- vacuum matrix element 



{0\A;{Q)\7T''{p))=t,Up^5^ 
extracted from the ratio of the two-point functions C^g and Cgg and 






(47) 



(48) 



in the large Euclidean time limit. 
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FIG. 7: Lattice QCD results for the primary ttAA coupling, CtAA- 
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FIG. 8: Lattice QCD results for the secondary ttAA coupling, HttAA- 
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The renormalized quark mass niq is determined from AWI, via two-point functions of the pseudoscalar density with 
either local (P^) or smeared {P^) quark fields. 



cEsit) = J2(^\T (p^(x,i)p3(o,o)) \n), 



(49) 



(and similarly for Cgg). The effective quark mass is defined by 



AWI/'^\ _ ™i- ^^ ^LSV) . l^SSV') 



'eff 



2 Zp Clsit) V Cisit) ' 



(50) 



and its plateau value yields mq. Note that Zp will be needed only if ones wants mq alone. Since Zp enters also 
Eq. (|13p . it cancels -as does Za since it comes with /,r- and therefore GttAA and Ht^aa are extracted directly from 
ratios of lattice three- and two-point functions without prior knowledge of either Za or Zp. We also note that the 
quark mass computed through ([50]) includes the effects of residual chiral symmetry breaking from the finite extent L5 
of the fifth dimension. These effects arc of the order of 17% for the DWF ensemble and 15% for the hybrid ensemble 
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K, or ami 


am,q af^/ZA Za 


Quenched Wilson fermions 


0.1554 
0.1558 
0.1562 


0.0403(4) 0.0611(14) 0.808(7) 
0.0307(4) 0.0587(16) 0.808(7) 
0.0213(4) 0.0563(17) 0.808(7) 


Hybrid or mixed action 


0.02 
0.01 


0.0324(4) 0.0648(8) 1.0994(4) 
0.0159(2) 0.0636(6) 1.0847(6) 


AT^ = 2 + 1 DWF 


0.004 


0.004665(3) 0.06575(12) 0.74521(2) 



TABLE II: The first column gives the hopping parameter k for Wilson fermions or the lattice mass of the domain wall fermion, 
the second the renormalized quark mass, the third the unrenormalized pion decay constant Jt^/Za in lattice units, and the 
fourth the axial current renormahzation constant Za- 

(also referred to as mixed scheme). Chiral symmetry breaking affects the PC AC relations and therefore the value of 
both strong couplings GttAA and H^^aa through Eq. ([13]). 

B. Testing Pion-Pole Dominance in the Axial and Pseudoscalar Matrix Element 

In this section we examine in detail the pion-pole dependence expected for the A form factors by performing fits 
to the results obtained. First, we test the validity of the Goldberger-Treiman relations of Eqs. (1^ [M)) by evaluating 
the ratios 

(51) 



MAQliq^) 

and 

MA/ii(g2) 



(52) 



These relations are expected to hold at low Q^. We show the results in Figs.lQlandfTOl The first ratio, given in Eq. (|5ip . 

carries moderate statistical error. It is consistent with unity for Q^ ^ 0.8 GeV^ for the quenched ensembles while it 
is underestimated at smaller Q^ values. This discrepancy at smaller Q"^ can be attributed to chiral effects on GttAA, 
which is expected to be more seriously affected by pion cloud effects than gi . The results using the hybrid dynamical 
ensemble, on the other hand, are consistently higher than unity for Q^ > 0.5 GeV^. The large statistical errors carried 
by these data make it difficult to draw definite conclusions. The behavior of this ratio is very similar to the behavior 
shown by the corresponding ratio for the nucleon GT relation as well as the nuclcon-to-A axial transition |21| . 

The second GT-type relation, given in Eq. (j52p . is statistically consistent with unity for the quenched results and 
Q^ > 0.8 GeV'^. The results from the dynamical ensembles are plagued by too large statistical noise to be able 
to meaningfully display them on the plot. We therefore have omitted these data from Fig. [TOl A very similar and 
consistent behavior with the first ratio is observed for the quenched data. We remind the reader that it is the first 
Goldberger-Treiman relation that is more significant for phenomenology, as it is this relation that connects the axial 
charge (from gi at Q^ = 0) to the G^aa coupling. 

To further probe the pion pole assumptions entering into our derivation of the GT relations we perform a set of 
fits to our form factor data. We have no a priori theoretical expectation for the functional form of 5i(g^), although 
typically a dipole form seems to accommodate well the nucleon axial form factor Ga as well as the leading axial N — A 
transition form factor C^. We note however that there seems to be a small dip in the gi at g^ = for the quenched 
ensembles. To accommodate this we fit the data to: 



(Q2 + m2) 



SiiQ') = ,1 ^,,, (53) 



The resulting fits are shown in Fig. [TT] The values for the fitted parameters are given in Table IIIII We note that 
the mass parameter mi determining the slope as Q^ — >■ is around 1 GeV, a scale typical for axial dipole masses 
controlling the dependence of nucleon Ga and the dominant axial TV — A C^ form factor. 
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FIG. 10: Ratio test of the Goldberger-Treiman Relation for HttAA- 



We consider the form 



a + bQ'^ 



{Q' 



(54) 



for 53 based on the pion pole dominance prediction given in Ea. (p5|) . The parameters o, b, and mi are fixed to the 
values arising from the fit of the gi data using the Ansatz given in Eq. (|53p . The fits are shown in Fig. [T^l The fitted 
parameters are given in Table IIVI Note that the value of m2 is considerably smaller compared to nii , as in fact is 
expected since this is detected from the presence of the pion-pole. This is especially verified by the quenched data 
where TO2 is close to the actual pion mass 111^^ of the ensemble. 
Pion pole dominance fixes completely the ratio gs/gi 
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(55) 
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FIG. 11: Fits to the data for the gi form-factor using the form given in Eq. 1531 



We form the ratio 53/(71 from our data and fit separately to a monopole form: 



(Q2+^2) 



(56) 



This fit is displayed in Fig.[T31 Using a ratio eliminates any need to know the theoretical form for gi{q^) alone. The 
fitted parameters c and 1712 are given in Table El The verification of the predicted form given in Eq. ([55)) is very good, 
with the pole mass 7712 consistent with the pion mass and the constant c reasonably close to 4Af^. 

The form factor hi is similar to 53 having a pion-pole dependence. We display the ratio hi/g^ in Fig. [14] for the 
quenched QCD ensembles. This ratio is notably constant over the whole Q^ range above 0.4 GeV^, with the constant 
~ 0.5. Based on this observation, wc use the Ansatz given in Eq. ([54]) also for hi. The fit is shown in Fig. [TSland 
the fitted parameters are given in Table [VD Again, TO2 is considerably smaller compared to mi, in accordance to the 
presence of a light (pion) mode. 

From Eq. ([2^ . the ratio h^/hi is completely fixed: 



hi 



AMI 



(57) 



We plot this ratio in Fig.[T6l Fitting the data to the monopole form of Eq. ([56| . we get parameters m2 and c within 
the range of the expected value (Eq. [57]) - see Table IVIII indicating that the subdominant form factor diverges with 
a double pion-pole- dependence. 

In Fig. [17] we present the fit of /13 to the Ansatz 



6Q2 



iQ^ + mjf 



(Q2 



(58) 
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FIG. 12: Fits to the data for the gg form-factor using the form given in Eq. 1541 



with a, b, and mi fixed to the values extracted from the fit of gi. The fitted parameters d and m2 are given in 
Table IVlIll in accordance to the hs/hi fit f Table IVlI)) . 

In the pseudoscalar sector, one expects a monopole dependence also for the ratio h/g. Fitting the data to the 
monopole form of Eq. (j56p , we get the parameters provided in Table IIXI Indeed, an agreement of the 7712 pole mass 
to the pion mass within the allowed by statistical noise regime is seen. 

The overall conclusion from the fits in this section is that all form factors satisfy qualitatively the pion-pole depen- 
dence predicted by PCAC. This is most clearly exemplified in the case of quenched QCD where the level of statistical 
noise allows such detailed analysis. In all cases the data fit these forms to good confidence levels, i.e., x^/dof < 1. 
Enhanced statistical noise for the dynamical ensembles limits the verification to the dominant form factors only, as 
the subdominant ones are beyond reach, but this still is a useful result as it shows the consistency between quenched 
and dynamical results. This corroborates other baryon studies that show small effects due to dynamical quark for 
pion masses larger than about 300 MeV. 
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FIG. 13: Monopole fits as given by Eq. [56]to tiie ratio gs/gi- 



niTr (GeV) a 



mi (GeV) x/dof 



quenclied Wilson fermions 

0.563 0.53(18) 2.15(31) 0.98(5) 0.82 

0.490 0.47(18) 2.08(33) 0.99(6) 1.07 

0.411 0.40(19) 1.98(38) 0.94(8) 1.48 

mixed action 

0.353 3.0(22.0) 2.4(1.9) 1.3(1.2) 0.44 

domain wall fermions 

0.297 0.19(24) 1.5(9) 0.82(18) 1.1 

TABLE III: Fit parameters for pi(Q^) using Eq. [53] 



m^ (GeV) 



ma (GeV) x7dof 



quenched Wilson fermions 
0.563 7.77(88) 0.54(11) 0.61 

0.490 7.45(89) 0.50(11) 0.55 

0.411 7.1(1.0) 0.44(15) 0.65 

mixed action 
0.353 10.7(4.6) 0.67(43) 0.49 

domain wall fermions 
0.297 10.0(5.3) 0.56(37) 1.26 



TABLE IV: Fit parameters for g-i{Q^) using Eq. [54] Parameters a, 6 and mi are fixed with the results of gi_ fits in Table Hill 
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FIG. 14: The ratio hi/gs as a function of Q^, with unity marked with a red line. 



m^ (GeV) ma (GeV) 



4Mi (GeV") x/dof 



quenched Wilson fermions 

0.563 0.523(64) 7.60(52) 8.64(18) 0.67 

0.490 0.477(63) 7.25(49) 8.12(18) 0.54 

0.411 0.396(75) 6.76(48) 7.64(21) 0.60 

mixed action 

0.61(18) 10.4(1.6) 9.40(33) 0.34 

domain wall fermions 
0.34(17) 8.6(1.3) 5.82(20) 0.66 



0.353 
0.297 
TABLE V: Fit parameters for g:i{Q'^)/gi{Q'^) using the monopole form of Eq. 1561 



m^ (GeV) 



7712 (GeV) x/dof 



quenched Wilson fermions 

0.563 3.04(48) 2 x 10"*(8 x lO"'^) 1.32 

0.490 3.32(84) 0.03(81) 1.37 

0.411 3.8(1.8) 0.1(1.4) 1.52 



TABLE VL Fit parameters for hi{Q^) using Eq. 1541 Parameters a, b and mi are fixed with the results of gi fits in Table Hill 



m,r (GeV) m2 (GeV) e 4M'i (GoV^) x7doT 

quenched Wilson fermions 
0.563 0.26(17) 7.6(1.1) 8.64(18) 0.21 

0.490 0.31(19) 7.4(1.1) 8.12(18) 0.38 

0.411 0.28(25) 6.7(1.1) 7.64(21) 0.58 



TABLE VII: Fit parameters for h3{Q'^)/hi{Q'^) using the monopole form of Eq. (56] 
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FIG. 15: Fits to the data for the hi form-factor using the form given in Eq. 1541 



m^ (GeV) d ma (GeV) xVdof 

quenched Wilson fermions 
0.563 22(6) 0.03(58) 0.46 

0.490 26(10) 0.25(28) 0.35 

0.411 28(15) 0.28(43) 0.33 



TABLE VIII: Fit parameters for h-iijQ^) using Eq. 1581 Parameters a, b and rrii are fixed with the results of gi fits in Table Hill 
The dynamical fermion data sets contain too much noise for the fits to be useful. 



m^ (GeV) ma (GeV) 



xVdof 



quenched Wilson fermions 
0.563 0.73(39) 4.2(1.6) 0.23 

0.490 0.42(43) 3.3(1.2) 0.22 

0.411 0.45(70) 3.9(2.0) 0.32 



TABLE IX: Fit parameters for h{Q'^)/g{Q^) using the monopole form of Eq. (56] 



20 



I silllii 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



30 

20 

10 



30 

=5" 20^ 
-^ 10 





30 

20 

10 





quenched m = 563 MeV 



ittittittuSi 



mi+ti+ti+Mmi+n^njg 



«|«|«i«i«i«m,^a»-i"-.TOi««^««,.ai^,^,„a^,„„^„„„ 




H \ \ \ \ \ h 



quenched m = 490 MeV 



WWHHaBWHHHH^HEH^H^. 



H \ \ Ph \ \ h 



quenched m =411 MeV 



WH^Hri:H^H^^^^,,,^^_^^_^^_^^ 



HHEH-tH4H. »jp. |-H-H-H-Hl [E gl 



J \ \ \ I I I I I I I I I \ \ \ \ \ L 



0.5 



1 

Q^ (GeV^) 



1.5 



FIG. 16: Monopole fits to the ratio hs/hi as described by Eq. 1561 
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FIG. 17: Fits to the data for the hs form-factor using the form given in Eq. 1581 



22 



m^ (GeV) gi(0) 

Quenched 

0.563(4) 0.589(10) 

0.490(4) 0.578(13) 

0.411(4) 0.571(18) 

Mixed action 
0.498(3) 0.573(23) 
0.353(2) 0.640(26) 

DWF 
0.297(5) 0.604(38) 



TABLE X: Numerical values for the dominant form factor gi (0) on each of the ensembles. Gaa 
ization. 



-3gi(0) with our normal- 



PHENOMENOLOGICAL COUPLINGS OF THE A AND COMBINED CHIRAL FIT 



Crucial parameters in heavy baryon chiral effective theories (HBxPT) with exphcit A degrees of freedom are the 
axial couplings of the nucleon, qa^ the axial iV — A transition coupling, CA^ and the axial charge of the A, ^aa. 
Assuming PCAC, these can be related via GT relations to the effective ttNN, ttNA and ttAA strong couplings: 



gA 



M. 



-g-rrNN 



CA 



N 






-gTvNA 



N 



.9AA 



A. 

Ma 



SttAA 



(59) 



We note that alternative notation and normalization factors exist in the literature in the definition of the effective 
strong couplings for ttA^A and ttAA. In addition, note that in such schemes Eqs. ((59)) are actually defining relations 
for the strong couplings. gA is very well known experimentally and a variety of lattice and theoretical calculations 
offer precise estimates, ca is much less-well determined, via the parity violating N-to-A amplitude which connects it 
to the dominant axial transition form factor C^{q^). gAA remains undetermined from experiment and is typically 
treated -as is also the case for ca- as a fit parameter to be determined from fits to experimental or lattice data. 
There have been several sum- rules calculations of the effective ttAA coupling [32M34l |. In Ref. [33 symmetry 

A_(_ and A_ form a chiral multiplet, lead to the conclusion that 



suggest that the 



arguments in a quartet scheme where N^, N* 

7rA±A± couplings (with like-charged As) are forbidden at tree-level. Quark-model arguments \li 

g-KAA = i4:/5)g^NN- 

Lattice calculations for the nucleon axial charge gA are available on a variety of ensembles and pion masses [26l | . 
In addition, results on the axial A^ — A transition from factor C^ [23 have been obtained on most of the ensembles 
used also in this work. We are therefore in position to perform a combined chiral fit using small scale expansion (SSE) 
within (HBxPT) [23, |3a, [33] • for gA, C^(q'^) and the A axial charge Gaa as functions of the pion mass 



The one-loop SSE expression for Cj has been worked by Procura [37 

rriTr is: 



The expression for C^{q ) as a function of 



C^ = ai + 02 ml + 03 g^ + loop5(m^), 



where the loop integral contribution is: 

CA 



loop5(m^) 



^ I ^ [ JffAA (407r ml + 101 A ml + 24A3) + ^ 5a Saa^^ A 



15552 7r2/^ I A L4""" ^ " " '2 

- 12A 4(162 m^ - 83A2) - 27 5^ (247rm^ -f- 75A?«2 - AQlS?) 

72 / 

+ -^ml-I^^(ml c\ - 28A2 c\ + \%g\ {ml - A^) 

50 



A 



A V 4 / \ m^ 

3 ml (900 c\ ~ ^.9aa - ^5aa gA + 81 5^ -f 648) 



50 
■8A2(-711ci-f-gi^ + 162.gi) 



In 



A 



(60) 



(61) 



Here A = A/a — Mn, A = 92.4 MeV, oi, 02, 03 are unknown parameters and A is a cutoff scale set to A = 1 GeV. 
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FIG. 18: Combined chiral fit: (a) Nucleon axial charge, qa, fitted to lattice data obtained with Nf = 2 twisted mass fermions 
(TMF) [2y|. Tlie pliysical value is shown by the asterisk; (filled circles: a=0.089 fm, L=2.1, filled squares: a=0.089 fm, 
L=2.8 fm, filled triangles: a=0.070 fm, L=2.2 fm, open square: a=0.056 fm, L=1.8 fm, star: a=0.056 fm, L=2.7 fm); (b) Real 
part of axial N to A transition coupling C5 (0) [22|; (open circles: a=0.124 fm, L=2.5 fm, filled square: a=0.124 fm, L=3.5 fm) 
and dynamical domain wall fermions (cross: a=0.114 fm, L=2.7 fm, open triangle: a=0.084 fm, 2.7 fm); (c) Real part of A 
axial charge Gaa ~ Sgi{0). 
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We use the SSE expression for the nuclcon axial charge is presented in Ref. [35 
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(62) 



(63) 



gr^ in the above expressions denotes the chiral hmit value of the axial charge, i.e. corresponds to gA in (|6T 
Finally, from Jiang and Tiburzi [27[ we obtain the chiral expansion for the axial charge of the A: 



GAA("^^) = ,9AA^A - 



(477/.) 



5AA'C(m., Ai) ( 1 + ^5aa 



+ C^ -SfAAA^(TO7r,-A,/i) -gAJ'(m.,-A, /x) 



The A field renormalization is 



1 



Ami 



1 
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Y^5AA'^("i^> m) + 2c^ J^(m., -A, 11) 



32^2/2 

and the loop integrals from ([331) evaluated at the scale n = 1 GeV: 

2 



£(m,^) 
/C(r?T,, A,/i) 



i^lo 



P-|a^),o,(^) 



-AVA2-TO2 1og f ^^ 1^ 

3 VA + yA2^ 



— — ;^ ( 77771 — V A^ — 7Tj2 fog I 



J^(m,A,/i) 



3 A 

l2 - 2A2 



A - VA2 - m2 
A + \/A2 - m2 



1 ( ^ 



2AVA2-TO2l0g(' 



VA2 - m2 



A + VA2 



(64) 



(65) 



(66) 



From the available lattice data on C^{q^ = 0;rn^), gA{m^) and G'aa(to2) we perform a simultaneous 7-parameter 
fit to expressions ([5n|) . (|62p and (|M)) fitting the unknown constants ai, a2, A, C^^^ as well as the common chiral 
couplings gA, ca and (7a A- Wc note that C^^^ is independent of ?n.; at a fixed value of A it can be fitted as a 
constant. 

The lattice nuclcon axial charge values g^(r7i2) are taken from twisted mass simulations [26| . Lattice values for 
the real part of the axial N — A couplings C^ (0) are taken from [22[ via a dipole extrapolation. The values of the 
real part of the axial charge of the A, G'AA('7i2 ) are related to the dominant axial form factor gi at zero momentum 
transfer via Gaa ~ ~3gi(0). For an additional lattice point to assist the fit we computed the zero-momentum gi 
values (only) on the 20'^ x 64 mixed-action ensemble with m. = 498 MeV. Values are provided in Table [X] 
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In Figure [TH] the combined fit is presented. The available lattiee data for all three observables vary mildly in the 
pion mass regime considered. gA remains underestimated with respect to the experimental value and the inclusion 
of C^ and Gaa into the SSE fit does not improve this systematically observed behavior. Strong chiral effects are 
expected at lighter pion mass values, especially below the A decay threshold, as is evident from the 1-loop trend of 
C^ and Gaa- 

VI. CONCLUSIONS 

A detailed study of the axial structure of the A(1232) has been presented, complementing recent and ongoing 
studies of the axial structure of the nucleon as well as the axial N~ to —A transition. The matrix element of the A 
state with the axial current has been parameterized via four Lorentz invariant form factors, gi, gs, hi and /13, and with 
two, denoted g and h, the pseudoscalar matrix element, generalizing the familiar nucleon axial structure. We detailed 
the lattice techniques required for the extraction of all six form factors for a complete q^— dependent evaluation via 
specially designed three-point functions. In fact, the calculation is optimized such that only two sequential propagators 
are needed for the numerical evaluation of the optimal correlators. PCAC constrains strongly the nucleon matrix 
elements of the axial-vector and pseudoscalar currents as is manifestly evident by the phenomenological validity of 
the Goldberger-Treiman relation. Lattice QCD provides a check of this relation, which is a result of chiral symmetry 
breaking present in the QCD Lagrangian, confirming that the g^-dependence of the axial and pseudoscalar form 
factors is in agreement with the PCAC predictions. Furthermore, the recent studies of the axial N— to —A transition 
have shown that PCAC constrains strongly also the transition from factors and measurements of the dominant form 
factor C5 and G^tNA provided a check of the non-diagonal GT relation. This work, examines extensions of similar 
relations for the A. The main result of the current work is that PCAC plays a major role also in the relation between 
the A matrix elements of the axial-vector and pseudoscalar currents, connecting the strength of the tt — A — A vertex 
to the A axial charge, Gaa- Actually, two independent pseudoscalar form factors, GttAA and /^ttAA are present and 
pion-pole dominance establishes relations among all six form factors. These predictions are qualitatively verified using 
results obtained in the quenched QCD study, which carries the smallest statistical noise. Results from two dynamical 
ensembles are consistent with these findings, albeit within large statistical errors. Having obtained an evaluation oi gA 
and C^ from previous studies and using the results of this work for Gaa on similar lattice ensembles, we performed 
a simultaneous chiral fit for all three utilizing one-loop chiral effective theory predictions in the SSE scheme which 
include a dynamical A field. The seven-parameter fit does not drive the prediction near the experimentally known 
gA value, and this is not surprising as it has become recently clear that the correct value of gA is not reproduced 
even with pion masses very close to the physical one. A careful isolation of excited state effects [40, |41| at pion 
mass of about 400 MeV failed to reveal excited state contamination in the lattice extraction of gA- Resolving such 
discrepancies is important for sharpening the predictive power of lattice QCD, which can yield phenomcnologically 
important quantities not accessible experimentally. Fully chiral 2-1-1 domain wall flavour simulations are available 
now below the 300 MeV pion mass utilized in this work, and this leaves open the perspective for further investigations 
in the future which will elaborate on the relations studied in this work and on the values of the major couplings that 
dominate the low energy hadron interactions. However, as shown here, the gauge noise is large and noise-reduction 
techniques will be needed in order to extract useful results using ensembles with close to physical pion masses. 
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Appendix A: Multipole Form Factors 



The axial vector transition between A states can be parameterized via a multipole expansion. This is most naturally 
performed on the Breit frame, where pf = —pi = (f/2. Let us denote the matrix element as 



(A(g72, Sf)\A ■ eA|A(-g72, s,)) = M(s/, s,;, A) . 
Generically four different transitions will occur parameterized via 

M(i,i,0) = L1+3L3 

M(||o) = 3L1-L3 

A/(i,-i,l) = -2E1-V6E3 

M(|^,l) = -V3E1+V2E3 



(Al) 



(A2) 



with Lj, Ej the longitudinal and electric multipole amplitudes of rank J. The polarization vector e\ has components 
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FIG. 19: Lattice results for the Ei multipole axial form-factor. 



We can relate the form factors gi, 173, hi and ^,3 to the multipole form factors Ei, £'3. Li and L^, which have physical 
relevance in the multipole expansion. 
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FIG. 20: Lattice results for the E^ multipole axial form-factor. 
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FIG. 21: Lattice results for the Li multipole axial form-factor. 



El = 
E3 = 

L3 = T 



^/2gl 2V2r(2ffi + fei(l+T)) 
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2T(2ffl+/ll(l+T)) 
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15\/rTT 
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Utilizing the above relations, we present results on the four multipole axial form- factors, Ei, E3, Li and L3, in 
Figures [H UHl [H and [H respectively. 

In the low momentum transfer limit, r ^ 0{(f' /M\) « 1 and from relations (|A5[IA7P we deduce that £'3, L3 << 1. 
On the other hand, £'1 and ii remain finite, as from relation (|A4p Ei ^ -\/2gi/3 and from (|A6|) ii ~ gi/3. Thus at 
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FIG. 22: Lattice results for the L3 multipole axial form-factor. 
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FIG. 23: Lattice results for the ratio of the Ei to V2L1 multipole axial form-factors for the quenched ensembles. 



the low momentum transfer limit we expect that Ei = V^Li + 0{(f). 

We test these predictions explicitly in Figure [23] where the ratio Ei/^/2Li is plotted. We observe a behavior 
consistent with a constant in the low energy (< O.SGeV ) regime although the numeric value of the constant is largely 
overestimated by the quenched lattice data. In addition, this constancy is in accordance to the pion-pole dependence 
of both El and Li which is evident from the quenched lattice data plotted in Figs. (|19p and (PT|) . Despite the large 
statistical uncertainties, E3 and L3 are consistent with small values at small momentum transfers. 
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We define Type I as 



Appendix B: Trace algebra for 3-point correlators 
1. Axial current correlator 

3 3 

i—1 cr.T—l 



(Bl) 



After evaluating the Dirac traces we find two distinct cases, fi = A and yu = 1, 2, 3. The kinematical frame is set to 
Pf = and p^ = -q. We note by ^ = (p^ + Af2^)i/2_ 
For /i = 4 we find 
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2. Pseudoscalar density correlator 

In a similar way wc evaluate the trace algebra for the pseudoscalar vertices. The index summation types are defined 
in the same way as in the axial case. Pseudoscalar Type I is 

3 3 

n|.s(<z) = E E ^-t^ [T'K.'{pf)Ol?,,KMp^)\ . (B9) 

z— 1 (7,r—l 

After the trace evaluation we find: 
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Pseudoscalar Type II is 
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after we evaluate the trace. 



Appendix C: Form Factor results 
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m^ = 563(4) MeV 



(GeV^) 
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Axial 
53 hi 
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0.562(27 
0.491(20 
0.459(19 
0.403(14 
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Pseudoscalar 
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2.12(14 
1.73(12 
1.35(13 
1.17(11 
1.05(14 
0.92(13 
0.73(10 
0.40(16 



46(89) 
46(17) 
17.4(8.5) 
10.3(5.5) 
10.1(2.9) 
8.6(1.9) 
5.5(1.3) 
4.21(92) 
2.34(90) 
2.08(69) 
1.95(72) 
1.52(57) 
1.18(43) 
0.48(60) 



TABLE XI: Delta form factors from quenched Wilson fermions. 



m^ = 490(4) MeV 



Q2 (GeV^ 



0.000000 
0.1728690 
0.3389351 
0.4989433 
0.6535119 
0.8031605 
0.9483310 
1.226705 
1.360524 
1.491113 
1.618694 
1.743467 
1.865609 
1.985279 
2.217769 



51 



Axial 
53 hi 



fts 



0.578(13 
0.695(72 
0.550(36 
0.474(28 
0.467(23 
0.408(17 
0.363(16 
0.326(17 
0.278(16 
0.252(17 
0.226(15 
0.202(20 
0.216(19 
0.173(16 
0.110(24 



14.0(2.9) 

7.03(78) 

4.52(45) 

3.99(31) 

2.88(18) 

2.18(14) 

1.58(12) 

1.253(95) 

1.14(10) 

0.935(78) 

0.71(10) 

0.751(87) 

0.580(71) 

0.41(12) 



20.0(8.3) 
5.1(2.0) 
2.4(1.1) 
2.45(60) 
1.53(33) 
1.14(25) 
1.02(18) 
0.67(14) 
0.42(13) 
0.40(11) 
0.54(14) 
0.501(87) 
0.348(73) 
-0.09(13) 



390(330) 

97(41) 
33(17) 
26.8(9.0) 
13.1(3.6) 
8.1(2.4) 
5.1(1.4) 
3.2(1.0) 
2.59(93) 
2.02(63) 
2.17(77) 
1.74(48) 
1.25(36) 
-0.02(58) 



Pseudoscalar 
5 h 



10.85(62) 

7.69(27 

5.78(23 

4.23(22 

3.54(17 

2.98(16 

2.21(17 

1.75(14 

1.30(14 

1.15(13 

1.02(16 

0.99(15 

0.73(12 

0.31(18 



120(130) 
59(23) 
26(11) 
10.3(6.8) 
11.6(3.3) 
8.9(2.3) 
5.9(1.5) 
4.1(1.1) 
2.0(1.1) 
2.01(82) 
1.50(89) 
1.79(70) 
1.10(50) 
0.21(72) 



TABLE XIL Delta form factors from quenched Wilson fermions. 



(GeV^ 



51 



Axial 
53 hi 



hz 



Pseudoscalar 
_J h 



m-n = 411(4) MeV 



0.000000 

0.1726468 

0.3381149 

0.4972318 

0.6506769 

0.7990168 

0.9427295 

1.217848 

1.349909 

1.478675 

1.604379 

1.727230 

1.847415 

1.965099 

2.193551 



0.571(18) 

0.77(11) 

0.546(52 

0.458(44 

0.495(41 

0.419(26 

0.361(23 

0.354(25 

0.276(22 

0.251(24 

0.224(20 

0.203(28 

0.230(24 

0.171(18 

9.100(30 



17.4(4.5) 
7.4(1.1) 
4.53(66) 
4.28(51) 
2.96(25) 
2.14(19) 
1.68(16) 
1.23(12) 
1.16(13) 
0.93(10) 
0.68(13) 
0.75(11) 
0.565(85) 
0.37(14) 



27(12) 
4.9(2.7) 
2.0(1.6) 
3.2(1.0) 
1.81(49) 
1.21(37) 
1.42(26) 
0.72(21) 
0.44(19) 
0.43(16) 
0.65(21) 
0.63(13) 
0.39(10) 
-0.18(19) 



560(490) 
102(61) 
33(24) 
32(15) 
13.4(5.2) 
7.5(3.4) 
6.4(2.0) 
3.4(1.4) 
3.1(1.3) 
2.26(87) 
2.5(1.1) 
1.87(64) 
1.32(47) 
8.70(83) 



13.1(1.1 
8.27(45 
6.17(36 
4.48(31 
3.73(23 
3.16(22 
2.41(25 
1.78(18 
1.23(18 
1.15(16 
0.96(21 
1.16(20 
0.73(15 
0.16(21 



290(230) 
59(40) 
27(18) 
15.4(9.8) 
14.4(4.4) 
11.2(3.2) 
7.5(2.3) 
4.4(1.4) 
1.8(1.5) 
2.4(1.1) 
1.0(1.3) 
2.9(1.0) 
1.13(68) 
-0.2(1.1) 



TABLE XIII: Delta form factors from quenched Wilson fermions. Results for hi and hs are plagued by statistical noise. 



33 



m^ = 353(4) MeV 



Q2 (GeV^ 



0.000000 
0.1240682 
0.2450231 
0.3630880 
0.4784607 
0.5913172 
0.7018154 
0.9162911 
1.020513 
1.122869 
1.223456 
1.322363 
1.419670 
1.515454 
1.702723 
1.794332 
1.884666 
1.973777 
2.061714 
2.148521 
2.234241 
2.402577 



91 



Axial 
53 hi 



0.640(26) 
0.62(39) 
0.51(17) 
0.41(12) 
0.28(10) 
0.385(69) 
0.349(60) 
0.368(57) 
0.283(47) 
0.222(50) 
0.234(45) 
0.212(66) 
0.262(46) 
0.237(48) 
9.586(67) 
0.148(44) 
0.151(39) 
0.190(66) 
0.141(42) 
0.173(56) 
0.06(31) 
0.14(39) 



15(26) 
9.8(6.0) 
6.0(3.1) 
1.9(2.1) 
4.0(1.1) 
3.42(78) 
3.22(60) 
2.07(43) 
1.39(43) 
1.43(36) 
1.31(51) 
1.50(31) 
1.32(30) 
0.41(38) 
0.70(24) 
0.65(20) 
0.82(32) 
0.66(20) 
0.77(26) 
0.4(1.6) 
0.5(1.5) 



8(74) 
-2(16) 
-4.4(7.7) 
-8.1(4.9) 
-0.9(2.5) 
-1.0(1.9) 
1.1(1.3) 
-0.6(1.0) 
-1.15(98) 
-0.50(77) 
-1.1(1.0) 
0.48(60) 
0.12(53) 
-1.15(83) 
-0.35(44) 
0.00(66) 
0.34(44) 
0.05(35) 
0.38(35) 
-0.1(2.6) 
0.7(2.3) 



-100(4500) 
120(540) 
-25(200) 
-160(104) 
-7(40) 
2(25) 
21(15) 
-0.3(9.3) 
-9.6(8.5) 
-3.2(6.3) 
-5.6(8.1) 
2.5(4.2) 
2.0(3.6) 
-6.4(4.8) 
-1.8(2.5) 
-0.3(2.0) 
1.3(2.3) 
0.5(1.8) 
2.1(1.7) 
0.0(11.0) 
2.5(9.2) 



Pseudoscalar 
g h 



14.4(4.8 
11.7(1.8 
9.1(1.2 
7.0(1.1 
5.88(67 
5.41(61 
3.46(56 
3.43(51 
2.22(52 
2.26(48 
2.34(77' 
1.70(50 
2.08(54 
1.25(99 
1.07(56 
0.90(48 
0.65(67 
0.37(52 
1.07(59 
-0.5(3.5 
1.0(3.4 



230(1800) 

440(320) 

260(140) 
87(83) 
85(37) 
87(27) 
29(17) 
35(12) 
17(11) 

13.3(8.5) 
18(13) 
7.0(7.1) 

12.3(6.3) 
14(11) 
8.4(5.2) 
3.1(4.0) 

-6.1(5.8) 
5.4(4.0) 

4.3(4.1) 
-5(32) 
2(18) 



TABLE XIV: Delta form factors from mixed-action fermions. 







Axial 




Pseudoscalar 




g' (GeV^) 


.91 33 hi 


h-i 


9 


h 




0.0 


0.604(38) — — 


— 


— 


— 




0.2060893 


0.08(45) -5.6(18.0) -44(45) 


-1500(2000) 


9.4(9.3) 


-1100(1800) 




0.4030337 


0.66(20) 9.3(4.3) 11(10) 


130(240) 


6.6(3.0) 


-300(280) 




0.5919523 


0.79(19) 10.3(2.9) 13(6) 


210(110) 


4.5(2.4) 


-82(140) 




0.7737530 


0.41(20) 3.5(2.3) 2.2(5.2) 


8(65) 


7.2(2.2) 


125(90) 




0.9491845 


0.21(12) 0.6(1.7) -1.0(2.5) 


-34(25) 


4.5(1.3) 


42(39) 




1.118873 


0.51(14) 3.9(1.1) 4.8(2.2) 


45(19) 


1.3(1.4) 


-20(32) 




1.443061 


0.47(21) 2.3(1.3) 4.3(2.6) 


21(17) 


3.2(1.9) 


27(30) 




1.598405 


0.31(14) 1.25(76) 2.1(1.6) 


4.7(8.8) 


0.1(1.2) 


-1(16) 




1.749719 


0.14(13) 0.91(69) 0.0(1.3) 


2.2(7.2) 


0.2(1.4) 


-12(16) 




1.897302 


-0.09(29) -0.3(1.4) -2.2(3.2) 


-11(16) 


-4.6(4.6) 


-62(57) 


m^ = 297(5) MeV 


2.041417 


-0.03(33) 0.1(1.6) -0.4(2.9) 


3(14) 


0.2(3.1) 


6(33) 




2.182297 


0.17(23) 0.51(91) 0.8(1.8) 


4.6(7.7) 


2.8(2.7) 


9(18) 




2.320150 


-0.06(14) -0.41(57) -0.7(1.1) 


-3.6(4.4) 


0.4(1.4) 


2(11) 



TABLE XV: Delta form factors from Domain- Wall fermions. 



